Fine Structure in the Off-Resonance Conductance of Small Coulomb-Blockade 

Systems 
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We show how a fine, multiple-peak structure can arise in the off-resonance, zero-bias conductance 
of Coulomb blockade systems. In order to understand how this effect comes about one must aban- 
don the orthodox, mean-field understanding of the Coulomb blockade phenomenon and consider 
quantum fluctuations in the occupation of the single-particle electronic levels. We illustrate such 
an effect with a spinless Anderson-like model for multi-level systems and an equation-of-motion 
method for calculating Green's functions that combines two simple decoupling schemes. 

PACS numbers: 73.23.Hk 



Resonant tunneling through small, isolated, multi-level 
systems such as a quantum dot Q| causes a peak in the 
zero-bias conductance whenever a single-electron level 
coincides with the chemical potential (//) in the leads. 
Without Coulomb interaction a series of peaks with sepa- 
ration Ae (energy level spacing) would be observed when, 
for instance, the energy levels are lowered relative to /i 
by external gates. At each peak the number of electrons 
in the system increases by one, filling the correspond- 
ing single-electron state. Coulomb interactions, however, 
have drastic effects on the conductance: (i) On the one 
hand the current can be suppressed over a large range of 
gate voltages since incoming electrons may be strongly 
repelled by those already present in the system. This 
phenomenon is generally known as Coulomb blockade 
(CB). A mean- field type picture suffices to de- 
scribe such an effect: Each time an additional electron 
is added to a single-particle state in the system all the 
other energy levels are shifted with respect to their pre- 
vious values by an amount U which is related to the 
Coulomb repulsion. Therefore, if U S> Ae, a sparse se- 
ries of conductance peaks with separation U is expected 
instead of that with separation Ae. (ii) A more exotic 
phenomenon can take place when, at very low tempera- 
tures, the conductance between CB peaks (off-resonance 
conductance) is enhanced assisted by quantum fluctu- 
ations in the single-electron, degenerate (usually spin- 
degenerate) levels HHiJ. This effect is closely related 
to the Kondo effect which is well known in the litera- 
ture 0. Below the Kondo temperature |6), what de- 
termines which one of the above mentioned phenomena 
dominates the conductance properties of the system is 
the ratio Ae/r, where T is the coupling strength to the 
leads. For Ae/T 3> I, the mean- field picture is basically 
correct and CB physics dominates ||. For Ae/r <C 1, 
fluctuations can take over and so can Kondo-type physics 



Although the CB phenomenon has been experimen- 
tally well established in quantum dots (I) , to our knowl- 
edge, the latter one is yet to be observed in such sys- 
tems, in part, due to the extremely low temperatures 
required. In this work we show that the presence of quan- 
tum fluctuations in the occupation of the electronic lev- 
els can also be playing a significant role close to the limit 
where CB dominates. From our results we conclude that 
there is a possibility of observing an enhancement of the 
off-resonance current through small quantum dots which 
carries a direct fingerprint of the discrete single-particle 
levels. Specifically we address the problem on how the 
mean- field picture breaks down as Ae/r crosses over from 
3> I to > 1. We find that in the regime Ae/r > I quan- 
tum fluctuations in the single-particle levels create "dy- 
namical" channels which are available for transport. In 
this limit, these new channels give rise to a multiple-peak 
structure in the off-resonance conductance in addition to 
the smooth signature of virtual tunneling processes which 
is known as elastic cotunneling . 

We begin by considering a multi-level system con- 
nected to left and right leads described by the Hamil- 
tonian 

N L N L 
H = Y^^d\d % + E kc\ c k+ E U ij4 d i4 d i + 

i=l keR.L j>i=l 

E E Vi(kMdi+4c k \, (i) 
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where d\ (di) are the creation (annihilation) operators 
associated with the Nl single-particle levels in the sys- 
tem with energies ej, ct (c*;) are the ones for the lev- 
els in the left and right leads with energies Ek , and 
Vi(k) are the hopping matrix elements between them. 
The third term in the Hamiltonian takes care of the 
electronic correlations within the system. Such a term 
contains the necessary contributions (those of the type 
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density-density interaction) to study the most fundamen- 
tal aspects of transport through quantum dots. (Addi- 
tional terms might be added if one is interested in more 
detailed correlation effects ||, but this lies beyond the 
scope of this work.) Degeneracies like those due to the 
spin degree of freedom are not considered either (a high 
magnetic field may be implicitly assumed). 

The conductance g through the interacting system can 
be calculated with the formula B 



2e^ 
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dujf / FD (uj)^sm{tr['yij(oj)Gi j (oj)]}, (2) 



where f'pi){to) is the derivative of the Fermi-Dirac dis- 
tribution function, ~f(uj) is the hopping matrix defined 

by 7 »7 L M/[7» + 7 L HL 
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G(lu) is the retarded Green's function that must be cal- 
culated in equilibrium, and "tr" denotes the trace over 
the levels of the interacting region jl(| . 

Equation-of-motion (EOM) techniques for calculating 
Green's functions have been used in the past in the con- 
text of the Anderson model and, recently, also in the 
context of quantum dots. []5|L2 13 According to such a 
technique one may write 



uj^d])) = ({d^d}}) + (([d^H^.d])), 



(3) 



where Gij(uj) = ((di] d^)) . The higher-order Green's 
functions generated by the last term in Eq. ([|) must 
be approximated at some stage of the calculation to ob- 
tain a closed set of equations. The simplest way to 
do that is in the Hartree-Fock approximation (HFA) 
where the higher-order Green's functions generated in 
the first EOM cycle are decoupled in the following way: 
((dinj;d\)) « (rij)((di; d\)), with (rij) = (d^dj). In this 
approximation the Green's function takes the simple form 
G(iu) = [(Gf°(w))- 1 -53(w)]- 1 ) where G°(tu) is the diago- 
nal Hartree-Fock Green's function for the isolated system 
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and E(w) = £ (u;)+ £ £ '(u) is the total coupling self- 
energy. The Green's function projected on a level i de- 
pends on the occupation numbers (rij) of all the other lev- 
els which are usually calculated self-consistently. Now, in 
order to include dynamical processes (which will turn out 
to be relevant in certain limits) in the Green's function, 
one has to improve upon the HFA. The way to do that 
in the framework of the EOM method consists in gen- 
erating higher-order Green's functions from ((djnjjdj)), 
namely, ((c k d]di\ d\)) , ((c\didj\ d\)), {{d^nr, d\)) , and 
((cferij; d\)). Following Hewson and Zuckerman p[ , 
we neglect Green's functions like ((ckdldi; dj)) and 



((c\didj-,d\)) which contain unpaired operators. Physi- 
cally, this corresponds to considering only one-electron 
processes. By generating new Green's functions from 
the remaining terms and successive decouplings similar 
to those mentioned in the HFA, a closed set of equa- 
tions is obtained. We will refer to this approximation as 
Hewson-Zuckerman approximation (HZA) from now on. 
After some lengthy algebra, and excluding off-diagonal 
terms, |l5| ] we find the following expression for the re- 
tarded Green's function: 
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where additional terms containing products in U up to 
jjN l -i are p resell ^ ) are not shown here. As in 
the HFA, the occupancies (rij) and correlation functions 
(rij . . .ni) may be calculated self-consistently. 

The HZA, as well as many other self-consistent EOM 
approximations beyond the HFA, |^,|l^,^3) presents, how- 
ever, a serious drawback: It gives unphysical values for 
the occupation numbers. Figure ^ illustrates this short- 
coming in the simplest two-level case. As usual, |(| 
we take the total coupling self-energy independent of lu, 
equal for both levels, and purely imaginary, T,jj(u>) — 
—iT. We plot half the difference in the occupation num- 
bers of the two levels as a function of Ae/T when the 
chemical potential lies in between the singly- and doubly- 
occupied states (notice that this corresponds to plot- 
ting the magnetization for the usual symmetric Anderson 
model djll]]). It is well known from other methods that 
the fluctuations in both levels are asymptotically sup- 
pressed as one increases Ae/T and/or U/T, and that, 
eventually, the HFA occupation numbers are basically 
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correct. ||,|6|,[l6|,|l7j For instance, second-order perturba- 
tion theories in U on top of the Hartree-Fock solution 
give ((n 2 ) - (ni»/2 w 0.4 for Ae/r w 4 and J7 « 27rF, 
which is what we obtain in the HFA. By contrast, the self- 
consistent HZA gives a much smaller value. Moreover, 
this value, instead of increasing with U /Y, decreases. One 
way to get around this difficulty and to obtain realistic re- 
sults out of the dynamical expression (^) is to avoid the 
fully self-consistent procedure. This can be effectively 
done by using the HFA static results for the occupa- 
tion numbers and multiple-particle correlation functions 
in the expression (^) for the Green's function (we will 
hereafter call this approximation HFA-HZ A) . In this way, 
as we approach the isolated-system limit, r — > 0, the cor- 
rect occupancies are guaranteed while keeping open the 
possibility for fluctuations at finite hopping. We expect 
this approximation to give reliable qualitative results in 
the limit Ae/r > 1. It cannot be valid, however, for 
Ae/r < 1 for several reasons: (i) The HFA breaks spon- 
taneously the local symmetry when U/Y > 1 (see Fig. 
|l|) ; (ii) off-diagonal elements have been ignored, and (iii) 
in order to reliably account for the strong fluctuation ef- 
fects, which occur in that limit, one should have gone 
beyond the HZA in the EOM method jLlJQ. 

We now calculate the conductance of a five-level sys- 
tem where an increasing coupling of the levels to the 
leads, = —i(Y x j), simulates a realistic situation 

for quantum dots. All the interaction terms C/y are set 
to one. Figure [2] shows g vs. p in two different limits. 
When Ae/r > 1 [Fig. |(a)] both HFA and HFA-HZA 
give similar results: There is a peak in the conductance 
whenever fi « E(N) - E(N - 1) = p do t(N) where E(N) 
is the ground state total energy of N particles in the 
dot. The single-particle states are successively filled and 
remain fully occupied as p moves up between the renor- 
malized single- particle levels. As expected, due to the 
derivative of the Fermi-Dirac distribution function in Eq. 
||, the heights and widths of the peaks are proportional 
to 1/T and T, respectively. In contrast, the level occu- 
pancies obtained from the HZA when off resonance [away 
from the charge-degeneracy points where p « Pdot(N)] 
are non-integer numbers and the total charge in the dot is 
not perfectly quantized. (Detailed analysis will be given 
elsewhere. Jl8[) 

As r increases [Fig. ||(b)] the occupation of the single- 
particle levels is no longer either strictly zero or one 
as a function of p. The mean-field picture rendered 
by the HFA is only approximately valid and quantum- 
mechanical fluctuations mediated by the interaction play 
their role now. In fact, within the HFA-HZA, new peaks 
appear in the off-resonance conductance as a consequence 
of these fluctuations. The presence of the peak labeled 
(1) at p = 1.0 can be understood in the following way: 
There is a finite probability for the electron to be in any 
of the levels 1, 2, 3, 4, and 5. It spends most of its 
time in level 1, but if it happens to be in levels 2, 3, 



4, or 5 a second, external electron can enter the sys- 
tem through level 1. Schematically these processes can 
be represented like | o • o . . .) — > | • • o . . .) — > | o • o . . 
oo». . .) — > |»o». . .) — > |oo». . .}, and so on, where empty 
and filled dots represent empty and occupied states at a 
given time, respectively. Although the cost in energy for 
these processes is the same, ej + e\ + Uij — e 3 = 1.0, 
their likelihood decreases with j. Feak (2) at p = 1.6 
can be understood in a similar way: There is a small 
but finite probability for the system to have the level 2 
empty even when p > pd ot (2). This is taken advantage 
of by an electron in the lead to sneak through via level 3. 
Schematically: | • o o . . .) — > | • o • . . .) — > | • o o . . .). The 
major peaks lie basically where the HFA predicts and all 
the other minor peaks can be associated with dynamical 
processes like those described above. [As can be seen, 
the self-consistent HZA overestimates the importance of 
these type of processes and gives rise to the spurious off- 
resonance structure seen Fig. 0(a).] 

It is worth mentioning that the tunneling processes 
described above seem to contribute to the off-resonance 
conductance in a way different from the usual elastic co- 
tunneling 0| which also anticipates a finite value of the 
off-resonance conductance. In such a theory, the contri- 
bution to the conductance comes from second- or higher- 
order (virtual) tunneling processes. Virtual transport is 
already included in the HFA where, in addition, diver- 
gencies close to the resonances are automatically taken 
care of. As can be seen in Fig. ||(b), the elastic cotunnel- 
ing is latent, for instance, in the slight asymmetry of the 
first major resonance. In agreement with previous work 
0, the Ae 2 /(Ae + Uf dependence of the off- resonance 
conductance can also be obtained in the HFA . In Fig. 
^|(b), however, we see that the off-resonance value in the 
HFA is always a lower limit of the one obtained in the 
HFA-HZA. 

Finally, we would like to stress the limits of validity 
of our results. At low enough temperatures, the off- 
resonance, fine structure relies on the existence of small 
quantum fluctuations in the occupation of the electronic 
levels. As we see from Fig. [l], such fluctuations are sup- 
pressed when U 3> Ae or Ae > T, but they survive if 
U > Ae and Ae > Y. Whereas the latter condition can 
be easily obtained by tuning gate voltages, the former one 
can only be generically found in relatively small quantum 
dots (diameter typically < 0.1pm). (Such quantum dots 
exist at present but their transport properties have not 
been fully analyzed to date.) When Y w Ae one expects, 
not only the fine structure to disappear, but the CB phe- 
nomenon altogether |nj (at least for a system with an 
infinite number of levels). This can be understood very 
easily. CB disappears when the conductance of the in- 
sulating barriers, which keep the dot isolated from the 
leads, becomes of the order e 2 /h. This is equivalent to 
saying that Ypdot ~ 1, where pdot is the density of states 
of the dot. Since, to a first approximation, p do t ~ 1/Ae 
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we obtain T/Ae w 1 as condition for the disappearance 
of the CB peaks. This trend is also confirmed by our cal- 
culations in the HFA-HZA for five levels as can be seen in 
Fig. ||. As the broadening of the single-particle levels ap- 
proaches Ae, the charge quantization is lost progressively 
and, consequently, the distinctive CB peaks change into 
a smooth oscillation as a function of fi. However, as we 
discussed before, a quantitative analysis of this regime is 
beyond the scope of our present work. 

In conclusion, an EOM that combines two different de- 
coupling schemes has been used to calculate the conduc- 
tance through a general multi-level system with strong 
interactions. A novel multiple-peak structure is obtained 
in certain limits. Whereas the main peaks are easily 
understood in terms of the orthodox Coulomb Blockade 
theory, quantum fluctuations must be invoked to explain 
the smaller peaks. The structure created by these small 
peaks in the off-resonance conductance lies on top of the 
contribution coming from virtual elastic tunneling. 
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FIG. 1. Difference in the occupation numbers of the two 
levels in the symmetric case as a function of Ae/r. The HFA 
and the HZA are shown for two values of U. The HFA gives 
the correct limit for large values of Ae/r when compared with, 
e.g., second order perturbation theory in U (see text). In the 
HZA, however, the difference is too small and decreases with 
U instead of increasing. 
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FIG. 2. Conductance in the HFA-HZA (thick solid line), in 
the HFA (thin solid line), and in the HZ A (dotted line) as a 
function of /j, in a five-level system up to N — 3. We have 
set ei = 0. The other parameters are: Ae — 0.3, Uij = 1, 
(a) kT = 0.02, T = 0.002, (b) kT = 0.002, F = 0.02 (all 
the magnitudes are in units of U). The smaller peaks in (b) 
correspond to dynamical channels opened by the fluctuating 
occupation numbers of the individual single-particle levels. 




FIG. 3. Conductance in the HFA-HZA and total charge 
as a function of fj, in a five-level system for different values 
of r (no dependence with the single-particle level has been 
considered now). kT = 0.005 and all the other parameters 
are as in Fig. @. As T increases, the fine structure starts 
disappearing as well as the overall CB effect. One can also 
see how the charge quantization is lost progressively. 



5 



